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We report the construction of a three-dimensional template bank for the search for gravitational waves from 
inspiralling binaries consisting of spinning compact objects. The parameter space consists of two dimensions 
describing the mass parameters and one "reduced-spin parameter", which describes the secular (non-precessing) 
spin effects in the waveform. The template placement is based on an efficient stochastic algorithm and makes use of 
the semi-analytical computation of a metric in the parameter space. We demonstrate that such a template family is 
effectual (effective fitting factor ~ 0.93-0.99) for the detection of generic spinning "low-mass" (mi +1112 < 12 M ) 
binaries in the advanced detector era, over the entire parameter space of interest (including binary neutron stars, 
binary black holes, and black hole-neutron star binaries). This provides a powerful and practical method for 
searching for gravitational waves from generic spinning low-mass compact binaries. The expected improvement 
in the average detection volume (at a fixed signal-to-noise-ratio threshold) of this spinning search, as compared 
to a non-spinning search is ~ 19 - 58%. 



I. INTRODUCTION 



The quest for the first direct detection of gravitational waves 
(GWs) is entering a new era. The first generation interfer- 
ometric GW detectors (LIGO, Virgo, GEO 600 and TAMA) 
have been decommissioned and are being upgraded to their ad- 
vanced configurations. Although a detection of GWs is yet to 
be made, the non-detection is consistent with our expectation of 
the event rates of the astrophysical phenomena producing GWs 
of detectable strength (see, e.g. [1]). The second-generation 
detectors are expected to be operational in a few years and to 
reach their design sensitivities within this decade. With roughly 
an order of magnitude improved strain sensitivity as compared 
to their first generation counterparts (and hence three orders 
of magnitude increase in the volume of the local universe ac- 
cessible to GW observations), the second-generation detectors 
are expected to make the first detections, opening up a new 
observational window to the Universe. 

Among the most promising sources for the first detection 
of GWs are the coalescence of astrophysical compact binaries 
consisting of neutron stars and/or black holes. The rationale 
for this expectation is at least threefold: 1 ) Compact binaries 
are efficient sources of GWs. During the coalescence process 
about 1-15% of the rest mass of the binary will be radiated as 
GWs, which means that such sources can be observed up to 
very large distances. 2) Radio observations of binary pulsars 
provide strong observational evidence supporting the existence 
of at least one class of such sources, i.e., binaries of two neutron 
stars. 3) These are remarkably "clean" sources — the expected 
GW signals can be accurately modeled and easily parametrized 
in terms of the component masses and spins. This means that 
the data analysis can benefit from quite sensitive detection 
techniques such as the matched filtering, which is the optimal 
filter to detect known signals buried in noisy data. 

Matched filtering involves cross correlating the data with the- 
oretical templates of the expected signals. Theoretical signal 
templates can be calculated by employing perturbative- or nu- 
merical techniques to solve the Einstein's equations of General 



Relativity 1 . During the early stages of the coalescence, called 
the adiabatic inspiral, the dynamics of the binary as well as the 
expected gravitational waveforms can be calculated using the 
post-Newtonian (PN) approximation to General Relativity [2]. 
For the case of "low -mass" binaries (total mass < 12 M ), the 
signal in the frequency band of the ground-based interferomet- 
ric detectors is dominated by the adiabatic inspiral. Hence 
it is sufficient to employ waveform templates modelling only 
the inspiral part of the coalescence, which is well described 
by the PN theory. On the other hand, for binaries with total 
mass > 12 M Q the "post-inspiral" stages (such as the merger 
and ringdown) also contribute to the signal observed by the 
ground-based interferometers [3, 4]. Thus, the templates have 
to model the complete inspiral, merger and ringdown. This 
requires inputs from the numerical-relativity simulations also, 
apart from perturbative General Relativity. 

Although the expected waveforms can be computed as a func- 
tion of the source parameters (such as the component masses 
and spins), the parameters of the signal that is buried in the 
data are generally not known. Thus, the data has to be cross- 
correlated with a "bank" of templates corresponding to different 
parameters. The discretization of the parameter space is gov- 
erned by two requirements: 1) the templates have to be placed 
with sufficient "resolution" in the parameter space such that the 
loss of signal-to-noise ratio (SNR) due to the mismatch of the 
closest template to a signal with arbitrary values of the parame- 
ters is minimal. 2) the computational cost of the search using 
the full template bank should be manageable. A geometrical 
formalism [5, 6] has been developed to lay down templates 
in the parameter space corresponding to a given value of ac- 
ceptable loss of SNR (or mismatch). This is based on defining 
a metric in the parameter space such that mismatch between 
two neighboring templates give the proper distance between 
them [5, 6]. 

While a closed-form expression of the template-space metric 
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can be computed for the case of binaries with negligible spins 
(or with spins aligned/anti-aligned with the orbital angular 
momentum), this is not possible for binaries with generic spins. 
This is due to the fact that, if the spins are misaligned with 
the orbital angular momentum of the binary, the spin-orbit 
and spin-spin interactions will cause the spins (and hence the 
orbit also) to precess. The resulting dynamics as well as the 
gravitational waveforms are rather complex, and the modelling 
requires solving a set of coupled ordinary differential equations. 
The template placement is further complicated by the large 
dimensionality of the parameter space (two mass parameters 
and six spin parameters). 

Almost all searches for GWs from coalescing compact bina- 
ries using the data of first generation interferometers assumed 
that the spins of the compact objects are negligible [7-12]. 
This was on the one hand motivated by the expectation that, 
for the case of first generation detectors, the non-spinning tem- 
plates are sufficient for the detection of spinning binaries over 
a significant fraction of the parameter space (loss of detection 
efficiency is acceptable). Another reason is the lack of a search 
strategy with improved efficiency (for a given false alarm prob- 
ability) as compared to the non-spinning search and that is 
computationally viable. 

The low-frequency sensitivity of the second-generation de- 
tectors is expected to be significantly better than that of the 
first generation detectors. For example, the Advanced LIGO 
detectors are expected to be sensitive to frequencies above 10- 
20 Hz, while the Initial LIGO detectors were only sensitive to 
frequencies above 40 Hz. Hence the second-generation detec- 
tors will be able to observe the inspiral from much larger orbital 
separations, and the observed GW signal will be significantly 
longer. As a result, neglecting the spin effects can cause a much 
larger dephasing of the template with the signal, and hence 
considerable loss of SNR [13, 14]. Proper consideration of the 
spin effects is essential in the advanced detector era. 

The fact that several different spin configurations are nearly 
degenerate has been recognized for quite some time. Different 
ideas for the construction of template banks for spinning bina- 
ries have been proposed in the past. All these proposals sought 
to reduce the effective dimensionality of the parameter space 
by making use of the near-degeneracies [15-25]; see Sec. I of 
[13] for a brief summary. Still, a computationally viable spin- 
ning search that is more efficient than a non-spinning search 
is yet to be demonstrated. Recently, it has been observed that 
binaries with moderate mass ratios (mxjm\ < 10) and generic 
spins can be detected using templates describing only the secu- 
lar (non-precessing) spin effects [13]. It was also demonstrated 
that these spin effects can be described by a single reduced- 
spin parameter in an approximate fashion. This opened up the 
possibility of performing searches for generic spinning binaries 
using a template family described by just three parameters (two 
parameters describing the masses and one describing the spins). 

In this paper, we extend the previous work presented in 
Ref. [13], which proposed a frequency-domain, spinning PN 
template family described by three parameters. Here we con- 
struct a computationally efficient 3-dimensional template bank 
for laying down templates in the parameter space. The template 
bank is based on an efficient stochastic template-placement al- 
gorithm and makes use of a fast, semi-analytical computation of 
the template-space metric for the waveform family mentioned 



above. Note that the template-placement algorithm (and the 
computer code) that is used for the construction of the tem- 
plate bank is quite generic that it can be used for constructing 
template banks for other waveform families as well. 

Finally, we note that other approaches are also being ex- 
plored towards developing efficient and feasible searches for 
GWs from spinning compact binaries. Other ideas include 
allowing the templates to take nonphysical values for the com- 
ponent masses, thus mimicking the effect of spins [25], sim- 
plifying the waveforms and the metric by neglecting the spin 
effects of the smaller body [17-22], reducing the effective di- 
mensionality of the parameter space by numerically identifying 
the principal components [26] etc. The number of templates in 
a bank can be further reduced by finding a near-orthonormal ba- 
sis of the template space by employing either the reduced-basis 
approach [27] or by employing singular-value decomposition 
of the template bank [28]. 

The rest of the paper is organized as follows: Sec I A presents 
a compact summary of the paper. Sec. II describes the reduced- 
spin PN waveform family describing the GWs from binaries 
with non-precessing spins. Computation of the template-space 
metric for this waveform family is described in Sec. Ill, where 
we also compare the analytical computation of the metric with 
exact numerical computations. Sec. IV discusses the construc- 
tion of the template bank using a stochastic placement algo- 
rithm and the template-space metric. The effectualness of 
the template bank in detecting generic spinning binaries is 
demonstrated in Sec. V, while Sec. VI contain some conclud- 
ing remarks. We use geometrical units throughout the paper: 
G = c = l. 



A. Summary of the paper 

A quick summary of this paper is as follows: We present a 
three-dimensional template bank that is effectual in detecting 
GWs from inspiralling compact binaries with generic spins. 
The parameter space consists of two dimensions describing 
the mass parameters and one dimension describing a "reduced- 
spin parameter" that describes secular effects of spin. The 
fitting factor of the reduced-spin template bank is 0.93-0.99 
for binaries with total mass < 12 M Q (see Fig. 4), while the 
corresponding fitting factor of a non-spinning template bank 
is 0.84-0.89. A search using the reduced-spin template bank 
is expected to bring about 19-58% increase in the detection 
volume at a fixed SNR threshold (see Fig. 5). 

The high effectualness of the reduced-spin template bank 
(which does not seek to model the modulational effects of 
precession) is due to the fact that, for the case of comparable- 
mass binaries (milm\ < 10) the total angular momentum of 
the binary is dominated by the orbital angular momentum, and 
hence the modulational effects of spin precession on the orbit, 
and hence on the observed signal is small. This effect is further 
enhanced by the intrinsic selection bias towards binaries that 
are nearly "face-on" with the detector (where the modulational 
effects of precession are the weakest while the signal strength 
is the strongest) as opposed to binaries that are nearly "edge- 
on" (where the modulational effects are the strongest while the 
signal is the weakest). 

The template placement is based on an efficient stochastic 



3 



algorithm, and makes use of the semi-analytical computation 
of a metric in the parameter space, which significantly reduces 
the computational cost. For the range of parameters that we 
have considered (see Table 1) the reduced-spin template bank 
results in a factor of ~ 10 increase in the number of templates 
(as compared to a corresponding non-spinning bank). This 
number is indicative of the expected increase in the computa- 
tional cost of the spinning search. We also emphasize that the 
template placement algorithm and the computer code (called 
SBank) is quite general and can be applied to arbitrary template 
waveforms with arbitrary dimensions. 



II. REDUCED-SPIN WAVEFORM TEMPLATES FOR 
INSPIRALLING COMPACT BINARIES 

During the early stages of the coalescence of the compact 
binaries, there is a clean separation of time scales. The orbital 
time scale (the orbital period) is much shorter than the pre- 
cession time scale (the time scale over which the spins/orbit 
precess around the total angular momentum axis), which is 
much shorter than the inspiral time scale (the time scale over 
which the orbital separation decreases), i.e., 



^oibit *^ ^precession *^ ^inspiral- 



(2.1) 



This clean separation of time scales considerably simplifies the 
equations to be solved for computing the expected GW signals. 
The gravitational waveforms can be computed by solving the 
following set of coupled ordinary differential equations [29]: 



dE(v) 

dt 
dSj 

dt 

c/Ln 

dt 



T(v) 



(energy balance) 



fi, xS;, i = 1,2, (spin precession) 
-1 d 



IILII dt 



(S1+S2). (orbit precession) (2.2) 



The first equation above is the energy-balance argument, which 
relates the PN expansions of the specific binding energy E(v) 
of the orbit with the GW luminosity T(v). The PN expansion 
parameter u is related to the orbital frequency coby v = (mcj) 1 ^, 
where m = rti\ + mi is the total mass of the binary. The second 
equation represents the precession of the spin vectors Si and 
S2, where the magnitudes of the vectors S~2 1 and $72 are the 
spin precession frequencies. In the leading order, and $^2 
are parallel to the Newtonian angular momentum vector Ln- 
The third equation represents the precession of the orbital plane 
(described by the unit vector Ln along the Newtonian angular 
momentum), which is derived from the fact that the total an- 
gular momentum is conserved over the precession time scale. 
||L|| is the magnitude of the orbital angular momentum L. Due 
to the precession of the orbital plane, and since the GWs are 
predominately beamed along the direction of the orbital angular 
momentum, the GW signal observed by a fixed detector will 
contain complicated amplitude and phase modulations. 

The essence of the reduced-spin templates proposed by [13] 
is the following: In the case of binaries with moderate mass 
ratios (miltni < 10), the total angular momentum of the binary 
is dominated by the orbital angular momentum. Hence the 



amount of orbital precession required to conserve the total an- 
gular momentum (by "compensating" for the spin precession) 
is rather small. As a result, the modulational effects of pre- 
cession on the waveform are small; the spin effects are nearly 
secular. Consequently, the phase evolution is very similar to 
that of a non-precessing binary (spins aligned/anti-aligned to 
the orbital angular momentum, so that dS,/dt — in Eq. 2.2) 
with a different value of the spins. This approximate mapping 
between precessing and non-precessing spins suggest that, we 
will be able to detect a lot of precessing binaries just employing 
non-precessing templates. This approximation (that the spins 
are non-precessing) has two advantages: 1) This enables us 
to compute an explicit, closed-form expression of the Fourier 
transform of the template using the stationary-phase approx- 
imation [30, 31]. 2) This enables us to describe the spin ef- 
fects using a single reduced-spin parameter, in an approximate 
way [13]. 

The "reduced-spin" waveform templates are defined in the 
frequency domain as [13]: 

h(f) = Cf- 1/6 exp {-i [¥(/) - */4]} , (2.3) 

where C is a constant that depends on the relative sky-position 
and orientation of the binary with respect to the detector, and 
/ is the Fourier frequency. Note that we have kept only the 
leading term in the frequency-domain amplitude. The phase of 
the GW signal is given by 
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where to is the time of arrival of the signal at the detector 
and 0o the corresponding phase, Vf is related to the Fourier 
frequency / by Vf = (nmf) 1 ^, m = ni\+ mi is the total mass 
and J] = mim^jm 1 is the symmetric mass ratio of the binary, 
and je is the Euler gamma. The spin effects in the waveform 
are known only up to 2.5PN order (u 5 ), and are described by 
the following parameters: 
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(2.5) 

where x is called the reduced-spin parameter, which is related 
to the individual spins of the binary by 



X = Xs + fyn 
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where 5 = (ni\ - m-^jm is the asymmetric mass ratio, and the 
symmetric- and antisymmetric combinations of the spins: 
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III. COMPUTATION OF THE TEMPLATE-SPACE METRIC 



A. Overview of the metric formalism 



can be interpreted as a metric in the parameter space. Thus the 
mismatch between two neighboring templates has the interpre- 
tation of the proper distance in the parameter space [5]: 



1 -M = gijAA'AA 4 . 



(3.5) 



A convenient way of computing the template-space metric 
gij is by projecting the Fisher information matrix on to the 
subspace orthogonal to A extr [6]. The fisher information matrix 
(of normalized waveforms h(f)) is defined as: 



T ab = (d a h(f;X),d h h(f-X)) 



(3.6) 



where d a denotes a partial derivative with respect to the param- 
eter A a . The indices a and b take values from 1 to 5 (including 
both intrinsic and extrinsic parameters). The template space 
metric gy can be computed by projecting F a h on to the subspace 
orthogonal to A extr . i.e., 



g = r, - r 2 1 r 4 



(3.7) 



This Section provides a brief overview of the metric for- 
malism proposed by Owen [5] for laying down templates in 
the parameter space. The waveform template h(f) defined 
in Eqs. (2.3)-(2.4) is parametrized by a set of parameters 
A = {Ai ntr , Aextr) where Ai ntl - are the intrinsic parameters (such 
as m, 7] and^-) that are intrinsic to the binary, and Aj ntl are the 
extrinsic parameters (such as fn and <po. The match between 
two neighboring templates in the parameter space is defined by 

M(X, AA) = max AAmr (h(f; A), h(f; A + AA)) , (3.1) 

where the angular brackets denote the inner product inversely 
weighted by the one-sided power spectral density S /,(/) of the 
detector noise, called overlap: 



(a,b) 



a(f)b*(f) + a*(f)b(f) 
S h (f) 



df. 



(3.2) 



Note that, in Eq.(3.1), the overlap is maximized only over the 
extrinsic parameters. The lower frequency cutoff /o is typically 
determined by the detector noise, which rises steeply below 
/o due to the seismic noise. The upper frequency cutoff / cut is 
due to the PN approximation breaking down when the binary 
approaches close separations (typically taken as the frequency 
of the innermost stable circular orbit). Also, the "hats" de- 
note normalized waveforms: h(f) = h(f)/ ^{h(f),h(f)). Note 
that, since the computation of the match (Eq. 3.1) requires nor- 
malized templates, we can effectively set C — 1 in Eq. (2.3), 
and hence the extrinsic parameters describing the location and 
orientation of the binary do not appear in the problem. 

Since the match function has its maximum value of unity 
at AA = 0, we can Taylor-expand the match about AA = 0. 
Truncating the expansion at second order, we get 
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Above g is a matrix with elements gij, where i and j take values 
1 to 3 (intrinsic parameters only). Similarly T\,T2, T$, T 4 are 
the sub-matrices of the Fisher matrix: 



Ta = 



Tu 




Tn 


r 2 i 


r 23 


r 23 


r 3! 


r 32 


r 33 


r 4 i 


T42 


r 43 


r 51 


T53 


r 53 



r, ee 



r 3 = 



Tu 


r 15 


r 24 


r 25 


r 34 


r 35 


r 44 


r 45 


r 54 


r 55 



(3.8) 
(3.9) 



B. Choice of coordinate system 

It is convenient to compute the metric in terms of a new 
set of variables {6>o, f? 3 , # 3 s), which we call dimensionless chirp 
times [6, 32]. The advantage is that, in this coordinate system, 
the metric components are slowly varying over the parameter 
space. The chirp mass (M c = mrj 1 ! 5 ), the symmetric mass ratio 
77 and the reduced-spin parameter^ can be written in terms of 
#0, 03 and 3S as 
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AA=0 



The dimensionless parameters f? , 63 and f? 3 s are related to the 
familiar chirp time [32] parameters to, t 3 and t 3 s by 

0o=27r/ O T O , 6» 3 = -2nf T 3 , 6 3S = 2nf r 3S , (3.11) 

where to is the Newtonian chirp time, t 3 and t 3 s are the spin- 
independent and spin-dependent terms of the 1.5PN chirp time, 
and /o is a reference frequency, such as the low-frequency 
cutoff of the detector sensitivity. 
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Using the relations given by Eq. (3.10), the phase of the 
reduced-spin waveforms given in Eq. (2.4) can be expressed in 
terms of O , 03 and O^s as 
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(3.13) 



The advantage of using this coordinate system to describe the 
waveform is that, at least the lowest order terms have near linear 
dependence on these parameters. This means that the derivative 
of the waveform with respect to these parameters, and hence 
the template-space metric is slowly varying over the parameter 



space (see, e.g., Figure 2). This makes the template placement 
problem considerably simpler. 



C. Template-space metric for reduced-spin templates 

Now we proceed to compute the template-space metric for 
the reduced-spin templates in the {0o, 03, 03s) coordinate system. 
For the efficient computation of the metric, it is useful to define 
the following noise moments: 
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(3.14) 



where f is the low-frequency cutoff of the detector, / cut is the 
upper frequency cutoff of the template (typically the frequency 
of the innermost stable circular orbit, and hence depends on the 
masses) and S/,(/) the one-sided power spectral density of the 
detector noise. The index p can take values between and 16, 
for a total of 17 values. 

These noise moments have to be computed only once during 
the construction of a template bank. Then the metric at different 
points in the parameter space can be computed using these noise 
moments. Using Eqs. (3.6), (2.3), (3.12) and (3.13), the Fisher 
matrix can be computed as: 

1 8 8 f 

Fa/, = - 2j 2j 1^0* d b i//( I k+e(f cut) 

k=o e=o ^ 

+ [d a if/\ d b i/f { + d a ifr k d b $] J k+t (/ cut ) 

+ d a iftd bt fi<K k+e (f cat )\ (3.15) 

where fa and 0r are given by Eq. (3.13). From the Fisher 
matrix T ab , the template space metric can be computed 
using Eq. (3.7). In this paper, we choose f cut as the frequency 
of innermost stable circular orbit of a test particle orbiting a 
Schwarzschild black hole: / cut = /sco = 6~ 3/2 (7rm)~ 1 . 



D. Comparison with numerical calculations 

Equation (3.3) provides a very good approximation of the 
match vV((A, AA) between two neighboring templates h(f; A) 
and h(f; A + AA). But this approximation becomes inaccurate 
for the case of two templates placed at large distances in the 
parameter space (due to the inaccuracy of the truncated Taylor 
expansion). In this section, we compare the approximate match 
function computed using the metric with numerically exact 
computation of the match, as given by Eq. (3.1). In the numeri- 
cal computation of the match, maximization of the overlap over 
the extrinsic parameters A extI = {to, 0o) is performed using the 
standard techniques — maximization over r is carried out by 
means of a Fast Fourier Transform, and the maximization over 
0o is carried out by employing two orthogonal templates [33]. 
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FIG. 1 : Comparison of the match ellipses computed from the semi-analytical calculation of the metric (black ellipses) with contours of the match 
function computed numerically (color contours). The component masses and reduced-spin parameter (mi,m 2 ,x) corresponding to point in the 
parameter space relative to which the match function is computed is shown on the top of each column (masses in units of M Q ). The different rows 
correspond to two-dimensional slices of these contours in the A6 - A8 3 plane (top row), A9 - A# 3S plane (middle row) and A6 3 - A8 is plane 
(bottom row). The solid black ellipses correspond to a match of 0.97 and dashed black ellipses correspond to a match of 0.99. 



Figure 1 shows the contours of the numerically computed 
match function (color-filled contours) along with the analyti- 
cally computed match contours employing the metric (black 
ellipses). For each value of {rn\,m%,x), the corresponding 
A = (#o,#3,#3s) is computed by inverting Eq. (3.10). The 
numerically computed match values of the template h(f; A) 
with neighboring templates h(f; A + AA) are reported by color- 
filled contours. Similarly, we compute the metric g,y(A) at 
the point A corresponding to (tn\,mi,x). Th e contour in the 
AA = (AOq, Aft), A#3s) parameter space corresponding to a 
match value of 0.97 can be found by inverting Eq. (3.3). This 
is an ellipsoid in this three dimensional parameter space. The 
black ellipses in Figure 1 are the two-dimensional slices of the 
ellipsoids. 

It can be seen from Figure 1 that the analytically computed 
match ellipses based on the metric agree quite well with the 
exact, numerical computation of the match function. The mi- 
nor disagreement between the two calculations is likely due 
to the fact that we are truncating the Taylor expansion of the 
match at quadratic order. But the approximate match ellipses 
are generally found to be smaller than the corresponding numer- 
ical exact match contours (this is consistent with the previous 
observations [6]). This essentially means that the template 
bank constructed using the metric will slightly over-cover the 
parameter space. 

Figure 1 also demonstrates the advantage of using the 
(8o, 03, #3s) coordinate system for laying down templates: Since 
the metric g is significantly flatter over the parameter space, 



the size and orientation of the match ellipses are very similar 
in all regions in the parameter space. In Figure 2, we plot 
yjgf, which is proportional to the local "density" of templates 
required to uniformly cover the parameter space. We see that 
yjgl varies by a factor of only as 10 over (6o, 03,@3s) space. In 
contrast, a metric computed in the (M c , rj,x) coordinate system 
has a variation of w 10 5 over the same parameter space (see 
Figure 2). The near-uniformity of the chirp-time coordinates 
is a desirable property for the speed of the template placement 
algorithm, which we discuss in the next section. 



IV. CONSTRUCTION OF THE TEMPLATE BANK 

We construct a template bank of the reduced spin waveforms 
using a stochastic placement method proposed in [34, 35]. The 
method begins with an initial set of intrinsic binary parameters 
Bo = { A; ntr , A^ ntr , . . . , A^ tr ), which is called the template bank 
seed and may be empty 2 . One constructs the full template 
bank by proposing randomly chosen values of the intrinsic 
parameters and checking whether the current template bank 
seed Bq already sufficiently covers the proposed template A'. 



2 Note that it is not necessary to provide discrete parameter values for the ex- 
trinsic parameters since the SNR of each template with the data is maximized 
over the extrinsic parameters using semi-analytical methods. 
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FIG. 2: The colored contours show the log of the square root of the determinant of the template space metric (log yjgj) over the range of template 
bank parameters chosen in Table I. The number density of templates in each region of the bank is roughly determined by yjgf. The top panel 
corresponds to the metric computed in the dimensionless chirp time coordinate system (8 , 83, 8 is ). The three subplots in the top row corresponds 
to three different slices in the 3-dimensional template bank, corresponding to^ = -0.3, 0, 0.3. The black dots correspond to the templates placed 
in each slice (of thickness Ax = 10~ 3 ) by the stochastic algorithm described in Section IV. Notice that the template density in different regions 
generally agree with the expectation from the metric. In the (8 , 83,83$) coordinate system, the density has a maximum variation of ~ 10 over the 
entire template bank. In contrast, we show in the bottom panels the same quantity computed in the coordinate system described by (M , n,x)- 
Notice that, in this coordinate system, the template density has considerable variation (~ 10 5 ) over the parameter space. This illustrates the 
advantage of using the (80, 83, 83s) coordinate system. 



We measure the coverage of the template bank via the fitting 
factor [36], which is computed by maximizing the overlap of 
the proposed template A' over the entire template bank Bq: 

FF = max M (A, A') , (4.1) 

\eB 

where M (A, A') represents the match of the seed template bank 
with the newly proposed template A' . The fitting factor thus 
gives the fraction of optimal SNR that can be obtained for the 
proposed template waveform A' using the existing template 
bank without including the newly proposed template. 

If the fitting factor for the proposed template is above a given 
minimum match threshold A1 m i n , then the proposed template 
is discarded to prevent over-coverage and we repeat the pro- 
cess with the same bank seed Bq. Otherwise, the proposed 
template is added to the bank and we repeat the process using 
B\ = Bq U {A'} as the new bank seed. The process continues 
until the number of consecutively discarded proposals exceeds 
a specified critical value fc max , indicating that the bank is statis- 
tically complete. 

The stochastic placement algorithm just described is straight- 
forward to implement and is completely generic, independent 
of many of the particular details of the waveform family one 
is using. For example, while it is useful to have a metric that 
is nearly constant over the parameter space, this is not a strict 
requirement for the stochastic placement algorithm, as opposed 



to lattice-based approaches [37]. We have implemented this 
algorithm with the generality of its use in mind. Thus, while in 
this paper we apply the stochastic bank placement code to the 
reduced-spin inspiral waveforms, the same code works with any 
waveform family, with or without an analytic approximation to 
the metric. We call the code SBank. SBank is implemented pri- 
marily in the Python programming language with speed-critical 
components in C. 

For practical applications, we labored to make SBank fast 
using some techniques suggested in related work [34, 35] and 
some novel. 

• We obtain algorithmic speedup by testing each proposal 
against not the whole bank, but only its neighborhood 
of templates, defined by some fractional difference in 
9q, the coordinate that is best fractionally measured; any 
template far away from the proposal in 9q cannot have a 
high match. 

• In the intermediate to late stages of bank construction, we 
will sift through thousands of highly matching proposals 
before finding one with small enough match to accept 
it into our bank. Short-circuiting this search early saves 
enormous computation. Thus we stop computing further 
matches immediately upon finding a match greater than 
the target minimum match. So that we find the high 
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Bank parameter 



Value 



Template waveform 
Noise PSD model 
Low-frequency cutoff: f 
Component mass: m\,rtii 
Total mass: m 
BNS reduced spin: x 
BH reduced spin: x 
Minimum match: A4 m in 
Max consecutive discards: k„ 



TaylorF2ReducedSpin 
aLIGOZeroDetHighPower 
20 Hz 
[1, 20] M 
[2, 21] M 
[-0.34, 0.34] 
[-0.9, 0.9] 
0.95 
10000 



TABLE I: Parameters used in generating the low-mass, reduced-spin 
template bank. The BNS region is where both mi and m 2 are less than 
2 M and the BH region is where either component mass is greater 
than 2M„. 



Bank 


Number of templates 


SBank reduced spin 


747977 


SBank non-spinning 


72029 


lalappsjtmpltbank non-spinning 


58484 



TABLE II: Comparison of bank size between our reduced spin bank 
produced by SBank with the parameters of Table I compared to the 
non-spinning SBank version and a matching lalappsjtmpltbank ver- 
sion. 



matches even sooner, we evaluate matches in the order 
of increasing 0q difference between the proposal and the 
bank seed. 

• Another technique we use is to draw proposals uniformly 
in (#0, #3, 6>3s) space. As the true template density is 
proportional to yfgf and g is slowly varying in these 
coordinates, this reduces the number of proposals thrown 
at already over-tested regions of parameter space and 
puts them in under-tested regions. 

• Finally, the availability of the metric gives an analytic ap- 
proximation to the mismatch, which significantly speeds 
up each iteration of this algorithm, but it is not strictly 
necessary. 

Using our stochastic template placement code SBank, we 
produced a template bank with a reduced-spin dimension de- 
signed to capture astrophysically plausible spins with the Ad- 
vanced LIGO zero-detuned, high-power configuration [38]. 
The template bank parameters are listed in Table I. The wave- 
form and noise model names refer to their designations within 
the LALSimulation software library [39]. 

We subjected this bank to a verification bank simula- 
tion, where we draw simulation waveforms from the same 
reduced-spin waveform family discussed in Sec. II (called Tay- 
lorF2ReducedSpin in LALSimulation) with parameter ranges 
given in Table I and we record the fitting factor of the bank 
towards each simulation waveform. The simulation parame- 
ters are chosen uniformly in (m\ , m.2,x) rather than (6q, 6*3, 63s). 
The results are shown in Figure 3. We see that in bulk, the 
algorithm has satisfied the minimum match criterion with only 
a very small leakage. 

Finally, in Table II we compare our reduced-spin template 
bank's size to non-spinning versions as generated by SBank 
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FIG. 3: Histogram of the achieved mismatch (1 - fitting factor) by the 
reduced-spin template bank in detecting reduced-spin injections. The 
plot is a demonstration of the achieved "coverage" of the bank. The 
vertical dashed line corresponds to a mismatch of 5%. Only ~ 1% of 
the ~ 5000 injections have mismatch > 5%. 



and lalappsjtmpltbank, the latter of which implements the 
hexagonal lattice-based algorithm used in the several most re- 
cent LIGO and Virgo inspiral searches [40]. In comparing 
the SBank non-spinning bank size to lalappsjtmpltbank non- 
spinning, we get a sense for the over-coverage of stochastic 
algorithms in two dimensions, which appears to be roughly 
23%. Comparing the SBank reduced-spin bank to the SBank 
non-spinning bank, we see that covering the spin space of 
interest requires ~ 10 times the number of templates, and ap- 
proximately the same factor in matched filtering computational 
cost. In the next section, we will see how much detection 
volume this extra computation allows us to access. 



EFFECTUALNESS OF THE TEMPLATE BANK IN 
DETECTING GENERIC SPINNING BINARIES 



TABLE III: Parameters used for the Monte-Carlo simulations of pre- 
cessing PN binaries. A binary component is deemed a neutron star 
(NS) if its mass is < 2 M Q and it is deemed a black hole (BH) if its 
mass is > 2 M G . 



Simulation parameter 



Value 



Waveform approximant 
BH spin magnitudes: \\xi\\ 
NS spin magnitudes: \\Xi\\ 
Cosine of spin orientations: L^'.x!" 
Location of the binary: (d L , 6, <f>) 
Cosine of the inclination: cos 1 
Polarization angle : iff 
Noise PSD model 
Low-frequency cutoff: f 



SpmTaylorTS 
uniform(0, 0.98) 
uniform(0, 0.4) 
uniform(-l, 1) 
uniform in volume 
uniform(0, 1) 
uniform(0, n) 
aLIGOZeroDetHighPower 
20 Hz 



In this section we study the effectualness [41] (a measure 
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0.89 0.89 0.87 



0.89 0.89 0.89 0.88 0.9 0.87 
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FIG. 4: The left plot shows the effective fitting factor FF c ff of the reduced-spin template bank at different regions in the component mass plane. 
Each filled circle corresponds to 5000 injections (with fixed component masses, corresponding to the center of each circle) of generic spinning 
binaries with parameters reported in Table III. The gray line shows the expected boundary above which the contribution of merger-ringdown 
becomes non-negligible (mismatch > 3%), and where the inspiral template bank needs to be replaced by an inspiral-merger-ringdown template 
bank. The black dashed lines correspond to the assumed boundary between the neutron-star- and black-hole mass. The right plot shows the same 
for a non-spinning template bank. 
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Injected mi (M ) 

FIG. 5: Average increase in the detection volume of a search employ- 
ing the reduced-spin template bank as compared to one employing a 
non-spinning template bank (corresponding to a fixed SNR threshold). 
The reduced-spin template bank is expected to bring about ~ 19 - 58% 
increase in the average detection volume. 



of the ability of a suboptimal template bank in detecting a 
family of target signals) of the reduced-spin template bank in 
detecting generic spinning binaries. We evaluate the effectual- 
ness by computing the fitting factor of the template bank with 
inspiral signals from generic spinning binaries spanning the 
astrophysically relevant parameter space. 

The fitting factor is the fraction of optimal SNR that can be 
obtained by using a sub-optimal template bank. This is com- 
puted by maximizing the overlap of the template bank (here the 
three-dimensional "reduced-spin" template bank /i templ (/; A)) 
with the "target" signals (here the generic spinning waveform 
family h taig (J; 9) ) over both the extrinsic A extl - = {to, <po} and 
intrinsic A; ntr = {m, rj,x} parameters describing the template 



waveforms: 

FF = max Acxlr>Am , r (h^\f; A), 1P*(f; 9)), (5.1) 

The calculation has to be repeated over different values of 
9, which describes the masses, spins and other parameters 
describing the relative location and orientation of the "target 
binary" with respect to the detector. 

The intrinsic luminosity of the target binary as well as the 
fitting factor of the templates depend not only on the masses and 
spins, but also on the parameters describing the location and 
orientation of the target binary. For example, the modulational 
effects of precession are the highest for binaries highly inclined 
with respect to the detector, while the intrinsic luminosity of 
such binaries is lower (as compared to binaries which are nearly 
"face on"). Thus, highly inclined binaries (which show the 
largest modulational effects of precession) are intrinsically less 
likely to be observed as compared to binaries that are face-on. 

In order to take into account such selection effects in eval- 
uating the effectualness of the template bank, we perform a 
Monte-Carlo simulation of generic spinning binaries and aver- 
age the fitting factor over the population. The waveforms are 
generated by solving the ordinary differential equations given 
by Eq. (2.2) in the "TaylorT5" approximation (see Sec. Ill of 
Ref . [ 1 3] for the full description) 3 . The target binaries (for a set 
of fixed values of component masses) are uniformly distributed 
in volume throughout the local universe. Spin magnitudes are 
distributed uniformly between zero and a maximum value (see 
Table III) and the spin angles are isotropically distributed. Co- 
sine of the angle i describing the relative orientation of the 
initial total angular momentum of the binary with respect to 



3 This particular approximant is chosen so as to disentangle the effect preces- 
sion from the effect of the difference between different PN approximants; 
see Appendix A for a discussion. 
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the line of sight is uniformly distributed in the interval (0, 1), 
while the polarization angle iff is uniformly distributed in (0, n). 
A summary of the parameters of the Monte-Carlo simulations 
is given in Table III. 

In order to evaluate the effectualness of the bank, we compute 
the effective fitting factor FF er r [17], in the following way: 



FF, 



eff 



3 ^ 
^bank 



1/3 



(5.2) 



where p = {h Ws , h mrs ) is the optimal SNR in detecting the tar- 
get binary, and pbank = P FF is the suboptimal SNR extracted 
by the template bank. The bars indicate ensemble averages over 
the full parameter space (while keeping the component masses 
fixed). The effective fitting factor FF e ff describes average de- 
tection range by a suboptimal template bank as a fraction of 
the detection range using an optimal template bank. The corre- 
sponding fractional detection volume (and hence the fractional 
event rates assuming that the binaries are uniformly distributed 
throughout the universe) is given by the cube of FF e( f . 

The estimated effective fitting factor FF e ff of the reduced-spin 
template bank is shown in the left panel of Figure 4. The figure 
suggests that the bank is effectual towards detecting generic 
spinning binaries over almost all the relevant regions in the 
"low-mass" parameter space {m\ + ni2 < 12 M Q ) 4 . Even in the 
worst case of the (4.5 M Q , 2.95 M Q ) binary, the effective fitting 
factor is ~ 0.93, and almost all regions in the parameter space 
the fitting factor is greater than 0.95 (note that the minimum 
match requirement M m i„ on the template bank was chosen to 
be 0.95). Note that the region above the gray line in the figure is 
the region where the contribution from the post-inspiral stages 
are expected to be significant, and the inspiral template bank 
needs to be replaced by an inspiral-merger-ringdown bank. 

The reduction in the fitting factor of the reduced-spin tem- 
plate bank in the high-mass and high-mass-ratio regimes is due 
to multiple reasons. The modulational effects of precession 
increase with increasing mass ratio, which are not modeled by 
our templates. There are additional factors causing the loss: 
The difference between different PN approximants become con- 
siderable at the high-mass-, high-mass-ratio regime (reflecting 
the lack of knowledge of the higher order spin-dependent PN 
terms), causing appreciable mismatch between the target wave- 
forms and the template waveforms even in regions where they 
should agree [42] (e.g., in the limit of non-precessing spins). 
Hence, it is likely that the fitting factor can be further improved 
by including the higher order PN terms, assuming that these 
higher order terms will reduce the difference between different 
PN approximants. 

The effective fitting factor of a non-spinning template bank 
(covering the same mass range) is shown in the right panel of 
Figure 4. The fitting factor of the non-spinning bank is 0.84- 
0.89 over the same parameter space. The average increase in the 



4 We conveniently define the "low-mass" range of the parameter space based 
on the previous studies using non-spinning inspiral waveforms, where it was 
shown that it is essential to include the effects of post-inspiral stages in the 
waveform for binaries with total mass > 12M Q [3,4]. 



detection volume provided by a search employing the reduced- 
spin template bank (as compared against the corresponding 
non-spinning template bank) is shown in Figure 5. The figure 
suggests that we can expect an increase of ~ 19-58% in the 
average detection volume at a fixed SNR threshold. Note that the 
real figure of merit of the improvement would be the increase in 
the detection volume for a fixed false-alarm rate. Calculation 
of this requires the calculation of the increase in the false-alarm 
rate due to the increased number of templates in the bank. We 
leave this as future work. 

Currently, all the observed neutron stars in binaries have spin 
periods > 22.7 ms [43], which correspond to spin magnitudes 
°f \\Xi\\ ^ 0.05. While this is not necessarily an upper limit 
on neutron-star spins, this could be indicative of the typical 
spins. We have repeated the simulations by restricting the spin 
range of neutron stars in the target binaries to the interval (0, 
0.05). This was found to make an appreciable difference only in 
the binary-neutron-star (m\ i < 2M ) region of the parameter 
space. In this region, the effective fitting factor of the non- 
spinning template bank was increased to 0.98. Thus, under 
this assumption, the non-spinning template bank appears to be 
adequate for the detection of GWs from binary neutron stars. 



VI. CONCLUSION AND OUTLOOK 

Developing an effectual and computationally viable search 
for inspiralling binaries of spinning compact objects has been 
a long-standing problem in GW data analysis. The problem 
is made difficult by the large dimensionality of the parameter 
space. In this paper, we have attempted one of the first, al- 
beit important, steps towards solving the problem: we have 
constructed a three-dimensional template bank that is effec- 
tual for the detection of a significant fraction of the generic 
spinning binaries in the "low-mass" parameter space. This 
development has been facilitated by a body of previous work: 
first, the realization that secular (non-precessing) spin effects 
are more important than the modulational effects for the case 
of comparable-mass binaries, which reduced the effective di- 
mensionality of the problem into three [13]. The computation 
of closed-form templates modelling GWs from non-precessing- 
spin binaries that are parametrized in terms of a "reduced-spin" 
parameter has made it possible to compute the template-space 
metric in a semi-analytic fashion. Secondly, the demonstration 
of computationally efficient stochastic -placement methods to 
place templates in the bank [34, 35] (which, as opposed to 
traditional lattice-based approaches, does not require the metric 
to be constant over the parameter space). 

We have demonstrated the expected effectualness of the tem- 
plate bank in the advanced detector era. For the spin distri- 
butions of the target binaries that we consider (see Table III), 
the effective fitting factor of the bank is in the range 0.93-0.99 
over the "low-mass" binary (m\ + m^ < 12 M Q ) parameter 
space. This is expected to bring about 19-58% increase in the 
detection volume as compared to a non-spinning template bank 
(for a fixed SNR threshold). The associated increase in the 
computational cost of the search would be roughly an order 
of magnitude. Note that further optimization of the template- 
placement algorithm and the parameter ranges is possible to 
reduce the computational cost. This first demonstration of a 
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template bank that is effectual (effective fitting factors > 0.93) 
over the entire parameter space of interest promises a pow- 
erful and feasible method for searching for generic spinning 
low-mass binaries (including binary neutron stars, binary black 
holes and black-hole neutron-star binaries) in the advanced 
detector era. 



Appendix A: Effectualness of the template banks against 
spinning waveforms generated in the TaylorT4 approximation 

In Sec. V, we studied the effectualness of the reduced-spin 
template bank in detecting a population of generic spinning 
binaries assuming that the target signals are given by the Tay- 
lorT5 PN approximant (see Sec. Ill of [13] for the full de- 
scription of this approximant). Note that this is one of the 
many approximations that can be used to compute PN wave- 
forms from inspiralling compact binaries, and these different 
approximants can produce somewhat different results (see [3] 
for an overview of different approximants). We would like 
to disentangle the loss of effectualness due to this effect from 
the loss due to the effect of precession. Thus, as the target 
waveform we need to use an approximant that is closest to 
the template in regions of parameter space where the target 
and template are expected to agree very well (e.g., in the limit 
of non-precessing spins). This is the motivation for choosing 
TaylorT5 approximant as the target waveform. 

We do not expect a priori one approximant to be closer to the 
signals given by nature than any other approximants. (This was 
further confirmed by comparisons of PN approximants with 
numerical-relativity simulations [44, 45]). Thus, in order to 
get a conservative estimate of the effectualness of the template 
banks, we compute their effective fitting factors towards sig- 
nals from generic spinning binaries computed in the TaylorT4 
approximation (Fig. 6). Note that the fitting factors are slightly 
lower than what we see in Fig 4. This difference arises from the 
fact that the waveforms computed using different approximants 
can be slightly different, reflecting the current uncertainty in the 
PN waveforms. It is likely that this uncertainty will decrease 
with the knowledge of higher PN terms (note that currently the 
spin-dependent terms are known only up to 2.5PN). 

There are additional factors contributing to the low effectu- 
alness of the template bank towards TaylorT4 injections. The 
particular implementation of the TaylorT4 waveforms that we 
use (called SpinTaylorT4 in LALSimulation) do not include 



the spin-dependent terms at 2.5PN order, while the reduced- 
spin templates (called TaylorF2ReducedSpin in LALSimula- 
tion) and TaylorT5 waveforms include these terms. Also, the 
TaylorT4 waveforms include the PN corrections to the time- 
domain amplitude, which are not included in the templates as 
well as TaylorT5 waveforms. These are only limitations of the 
particular waveform generation codes, and can be resolved by 
adding the appropriate PN terms. 

We argue that the main reason for the lower effectualness of 
the reduced-spin template bank towards TaylorT4 waveforms 
at high mass ratios is not the modulational effects of preces- 
sion; but the difference between PN approximants. In order 
to demonstrate this, we compute the effective fitting factor of 
the reduced-spin template bank towards TaylorT4 waveforms 
with non-precessing (aligned/anti-aligned) spins (Fig. 7). It 
can be seen that the effectualness of the template bank in this 
case is comparable to its effectualness towards generic spin- 
ning TaylorT4 waveforms (see left plot of Fig. 6). The loss 
of effectualness of the template bank in this case cannot be 
attributed to the effects of precession. These results greatly 
motivate the need of computing higher order spin terms in the 
PN approximation. 
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